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Visit the website of TIME 2012 
http://time2012.ut.ee/ 

The conference will be held on July 10-14, 2012 in Tartu, Estonia 


Fortunately four excellent experts have accepted our invitation to give an invited lecture: 

Bill Bauldry, USA (CAS in the Classroom: Yesterday, Today, and Tomorrow) 

Djordje Kadijevich, Serbia (Critical issues of effective CAS utilization) 

Walter Wegscheider, Austria (Development of pedagogical tools for teaching 

mathematics -from the PC Lab to Smartphone) 

Jose Luis Galan, Spain (From Vienna to Tartu: A 10 Years Tour using Derive in TIME) 

Additionally we are expecting an Estonian speaker, too 

Marina Lepp, Eno Tonisson, Josef Bohm (Co-Chairs of TIME 2012) 
(The cornflower = bluebottle is the national flower of Estonia, its picture is CAS-made!) 

Td like to invite you reading my address at the occasion of 20 years TIME-Conf- 
erences. You will find it on the TIME 2012 website, Josef 


Some interesting links (most of them provided by Michael de Villiers, SA) 

http://math.pppst.com/problemsolving.html 

http://www.qedcat.com/ 

http://www.lettredelapreuve.it/ 

http://www.scribd.com/doc/18935265/Geometrv-Unbound 

http://www.pdfbe.com/po/polvnome-book.pdf 

This is a collection of talks given at several “Didactic Days” of the Austrian Mathemat- 
ical Society (in German): 

http://www.oemg.ac.at/DK/Didaktikhefte/ 

The next two links are leading to (German) papers which deal with linear difference 
equations (David Sjostrand’s Proof): 


http://www.oemg.ac.at/DK/Didaktikhefte/2009 Band 42/VortragPauer.pdf 
http://homepage.univie.ac.at/david.wozabal/ws2004/diffgl.pdf 
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Liebe bUS-Mitglieder, 
nach vielen Problemen mit meinem Note- 
book und dem PC - Crash fur das erstere 
und Speicherprobleme mit dem zweiten - 
kann ich Ihnen endlich den neuen bNL#82 
vorlegen. 

Ziemlich gleichzeitig ist auch die neue 
Version des DNL#26 fertig geworden. 
Beide Newsletter waren sehr arbeitsin- 
tensiv. Wdhrend im „26er“ die Anpassung 
einiger DERIVE-Routinen an die Version 6 
manchmal einiges Kopfzerbrechen berei- 
tete, war es in der aktuellen Nummer 82 
die Obersetzungsarbeit fur die Beitrdge 
von D. Oertel und P. Trebisz. 

Ich schdtze sehr die genauen Erlduterun- 
gen von P. Trebisz und danke fur seine 
Bemuhungen. Mit Herrn Oertel pfiege ich 
nun seit Idngerer Zeit einen angeregten 
Briefwechsel und verfolge seine interes- 
santen unkonventionellen Ideen, die er uns 
im Lauf der ndchsten DNLs auch ndher 
br ingen wird 

Oavid Sjostrand verdanke ich eine weitere 
sehr intensive Beschdftigung mit dem 
„Proof' aus einem f riiheren DNL. 

Ich freue mich auch, hier wieder einen 
Beitrag von Michel Beaudin in seiner Mut- 
tersprache Franzdsisch verdffentlichen 
zu kdnnen. Michel hat - angeregt von Al- 
bert Rich's RUBI (Rule Based Integration) 
- seinen Unterricht hinsichtlich Integra- 
tion wesentlich verdndert. Er hat dariiber 
einen spannenden Vortrag im Rahmen der 
ACAll in Houston gehalten - den wir mdg- 
licher-weise auch hier im DNL werden 
verdffentlichen diirfen. 

Beachten Sie bitte die Hinweise auf TIME 
2012 sowie die website zum Download der 
TIME 2010 Proceedings (Seite 30). Jose 
Luis Salon, Sabriel Aguilera und Pedro 
Rodriguez haben auch hier einen perfek- 
ten Job hingelegt. IMuchas Sracias, ami- 
gos! 


Dear DUS Members, 

after a lot of problems with my laptop and 
my PC - the laptop crashed and the PC 
demonstrated its memory problems - I am 
glad to present the DNL#82. 

Pretty to the same time the revised 
version of DNL#26 was ready. Both 
newsletters needed a lot of work. In 
DNL#26 it was the adaptation of some 
DERIVE functions on DERIVE 6 and in 
DNL#82 much time was spent for trans- 
lating the contributions of D. Oertel and 
P. Trebisz. 

I appreciate very much the accurate 
explications of P. Trebisz. I am very 
grateful for his efforts making his work 
understandable for us. D. Oertel and I 
have been enjoying an interesting ex- 
change of mails since several months. I am 
following his unconventional ideas which 
he will present in the next DNLs. 

I owe David Sjostrand another very 
intensive activity in connection with “The 
Proof" from an earlier DNL. 

I am glad to present a contribution from 
Michel Beaudin in his mother tongue. 
Michel has - inspired by Albert Rich's 
RUBI (Rule Based Integration) - changed 
his teaching of integration. He gave an 
exciting lecture about his point of view at 
ACAll in Houston just recently. I hope 
that we can publish his paper in one of the 
next DNLs. 

Please notice the remarks wrt TIME 2012 
and the website for downloading the 
Proceedings of TIME 2010 (page 30). 
Jose Luis Salan, Gabriel Aguilera and 
Pedro Rodriguez did an excellent job. 
IMuchas Gracias, amigosi 



Download all DA^L-Derive- and Tl-files from 

http : / / WWW . austromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the 
Bulletin of the DERIVE & CAS- 77 User 
Group. It is published at least four times a 
year with a contents of 40 pages minimum. 
The goals of the DNL are to enable the 
exchange of experiences made with 
DERIVE, T7-CAS and other CAS as well 
to create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to 
reinforce the international touch of the 
DNL. It must be said, though, that non- 
English articles will be warmly welcomed 
nonetheless. Your contributions will be 
edited but not assessed. By submitting 
articles the author gives his consent for 
reprinting it in the DNL. The more contri- 
butions you will send, the more lively and 
richer in contents the DERIVE & CAS-T7 
Newsletter will be. 


Next issue: September 2011 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encription, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - whafs new? J. Bohm, AUT 
Truth Tables on the Tl, M. R. Phillips, USA 
Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 
Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ & Rob Gough, UK 

A Conics-Explorer, J. Bohm, AUT 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Treating Differential Equations (M. Beaudin, G. Piccard, Ch. Trottier), CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Statistics with Tl-Nspire, G. Herweyers, BEL 

Cesar Multiplication, G. SchodI, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 
Airplane Wings and Conformal Mapping, K. Korner, GER 
Rational Hooks, J. Lechner, AUT 
Conics with Tl-Nspire, W. Propper, GER 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, DXust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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This is the result of an inspiring exchange of emails between David 
and me and kept both of us busy for a while. Many thanks David for 
this exciting investigation, Josef 


The Proof, David Sjostrand 's Contribution 

I received an email from David Sjostrand (31 January 2011) 

Dear Josef, 

after only having had a short glance at your and Carl’s paper in DNL #80 1 couldn’t resist 
investigating some more sequences with Excel. (I shouldn’t have done this. 1 have so many 
things I must do.) 

It seems that the sequence P(n)=m*P(n-3) + (l-m)*P(n-2) is convergent for all m and all 
initial conditions. 

I attachjosef_carl.xls. 

Best regards, 

David 
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A 

B 

C 

D 1 

1 — ^ — 

F 

1 




'"i \ i 



2 

P = 

mP\ 


\-m\l 

Pn 7 


3 

71 

n 

—6 

^ 

71 — Z 


4 

34 

73 



in 

1-m 

5 

5 

3 



0,2 

0,8 

6 
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7 

10,8 

17 


1 

1 50,81818181 

1 -8,3636 

S 

81 

-28,2 



50,8181818' 

' -8,3636 

9 

28,64 

6,4 



50,8181818 

-8,3636 

10 

66,96 

-19,16 



50,8181818 

-8,3636 

11 

39,112 

-0,52 



50,8181818 

-8,3636 

12 

59,296 

-14,048 



50,8181818 

-8,3636 

13 

44,6816 

-4,248 


initial values 

50,8181818 

-8,3636 

14 

55,2592 

-11,342 



50,8181818 

-8,3636 

15 

47,6045 

-6,208 



50,8181818 

-8,3636 

16 

53,1437 

-9,9235 



50,8181818 

-8,3636 

17 

49,1354 

-7,2349 



50,8181818 

-8,3636 


Josef carl.xls 


My answer was (1 February): 

Dear David, 

thanks for this. It is an interesting idea. I reproduced your sequence(s) with DERIVE 
and wondered about the locus of the final (= limit) points with the first three points 
given and 0 < m < 1. 

I introduced a slider bar for m and plotted the 20*'^ points of the sequences with an 
increment of 0.05 for m - interesting. I will proceed investigating this. 
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Some days later another mail arrived (2 February): 


Dear Carl and Josef, 

thank you again for an interesting paper, J. Bohm & C. L. Leinbach: Using Rational Arithmetic... 
I think that I have a solution to the challenge on page 22. 1 attach some thoughts. 


^ =[«pA] 


Let < 


[«2,&2] 

[a3,&a] 


Pn={\^yn) = mPn-2 + ^Pn-3 » « > ^ 


I am almost sure that I can prove that 

m= — ^ hm/^ = —\a^ +2a^,b^ + 2bA 

using the theory of linear homogenous difference (recurrence) equations. 
It is easy to prove that lim = 0 or is not convergent if m^ — 

n^co 2 


Here is a happy, but very, very tired grandfather and two granddaughters, Josefine's daughters Cecilia 
and Emma. 



Best regards, 
David 


From: Yvonne Fogelander-Sjostrand <yvonne@ydsa.se> 

Datum: Thu, 3 Feb 2011 21:29:58 +0100 

Dear Josef and Carl, 

please find enclosed what I believe are proofs of the assertions I made in my last mail. 


Best regards, 
David 
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Px={aM 

Let ( 

— ycL^ 5 63 ] 

Pn={\’yn) = r^Pn-2 + ^Pn-i ^ ^ 

I am almost sure that I can prove that 

m = — => limP =— f«i +2aj + 2 a.,by +2&2 +2&,1 

using the theory of linear homogenous difference (recurrence) equations. 


It is easy to prove that lim = 0 or P„ is not convergent if m^ — 

n^cc 2 

A difference equation, + ax^_^ + bx^_^ + cx^_^ = 0 , has the general solution 


x^=Ar"+Br"+Cr", 


A, B and C are arbitrary constants, where r^,r2, are the roots of the characteristic equation 


+ ax^ + bx + c = 0 . 


(1) gives us 2 difference equations 


x^ - mx^_2 - tnx^_2 = 0 , x^= a^, x^=a2 and X3 = 


yn - W „-2 - ^yn -2 = 0 , Pi = , ^2 = ^2 ^nd ^3 = 63 ( 5 ) 

T 1 

1 . m = — 

2 

( 4 ) gives us x„ --^x„_2 = 0 , x^=a^, X2=a^ and X3 = (6) 

The characteristic equation x^-— x- — = 0 has the roots r, = , r, = -—-—and r, = 1 , 

2 2 ' 22 ^ 22 ' 

The general solution is x„ = Ar" + Br^" + C . Since |^i | < 1 and |^2 1 < ^ conclude that 

lim x„ = C . The initial condition x^=a^, x^= and X3 = gives us a possibility to find C. 


We find that C = same with the sequence , ( 2 ) is 

proved. 


II. m . 

2 

Suppose that the sequence x„ is convergent and that lim x„ = x . Then x-mx-mx = 0 which 

n^co 

gives us (1 - 2 m) X = 0 => x = 0 since . Thus ( 3 ) is proved. 
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From: nojo.boehm@pgv.at 

Sent: Friday, February 04, 2011 9:27 AM 


Dear David, 
dear Carl, 

Thanks for this elegant proof. If we had known this prior to TIME 2010 our talk would have 
lasted only a few minutes. 

I attach a DERIVE file dealing with your sequence from your last mail: 

P(n) = m*P(n-3) + (1-m)*P(n-2). 

I started investigating the locus of the limit points for 0 <= m <= 1 for any initial triple of 
points. 

Best regards to you all, 

Noor and Josef (still skiing in Hinterstoder = www.hiwu.at ) 


# 1 : 


Branch := Real 


ptsmCi'm'i m , n, pt) : = 

Pnog 

pt := I'm' 

k := 4 
#2 : Loop 

If k > n 
RETURN pt 

pt := APPEND(pt, [nfi-pt 4 .(k - 3) + (1 - m)-pt|.Ck - 2)]) 
k :+ 1 


These are the first 1 D points of your example sequence: 
34 73 

0 . 2 . 10 


#3 


#4: 


ptsni 


LL 100 

-36 J 

34 

73 

B 

3 

100 

-36 

10. S 

17 

SI 

-2S.2 

2S.64 

6.4 

66.96 

-15.16 

39.112 

-0.B2 

B9 . 296 

-14 . 04S 

. 44.6S16 

-4 . 248 


I am trying to plot an approximation of the locus of all "final" points (taking point #30 of each 
sequence) with m - which is k in the expression below - running from 0 to 1 with an 
increment of 0.02. 
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VECTOR 


#B: 


ptsm 


0 0 
B 3 
-2 -4 


, k, 50 


, k, 0, 1, 0,02 


30 


Inserting a slider for m with 0 < m < l (50 steps) I can observe the approximated limit point 
walking on the "locus". 


# 6 : 


ptsm 


0 0 
B 3 
-2 -4 


m , BO 


30 


Then I tried to insert sliders for the coordinates of the initial points. This didn’t work because 
of too long calculation times. I could realize this with GeoGebra. 



Another example 




VECTOR 


#7: 


ptsni 


0 0 
-B 3 

LL 2 10 J 


k. BO 


k. 0. 1. 0.02 


30 



3 
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We obtain a similar graph. 

My next question is: What is the locus of the points? 


I wrote that I could not realize moveable initial points with DERIVE. But I could do this with 
TI-NspireCAS (version 3). 


First of all I generated three 
points (the initial points of the 
sequence). The coordinates 
of the points must be defined 
as variables. After this I impl- 
emented a slider for m. 

You cannot see the red point 
- which is the 50*'^ element of 
the sequence - at this 
moment. 


I transferred the coordinates 
into the spreadsheet and 
created the first 50 points of 
the sequence. 


♦ 

A| 

B| 

C:i 




46 

47 

48 

49 

50 

50.8181 

-8.3636 


50,8182 

-8.36366 


50.8182 

-8.36362 


50,8182 

-8.36365 


50.8087 

-8.35726 


Cl 






Please compare the values with the Excel- and the DERIVE results. 

The coordinates of the 50*'^ point are made to variables, and these variables are used to plot 
the red point in the Graphs Application. 

Moving the slider one can observe the path of the “limit point” for a given set of initial points. 
You can also drag the points. 


From: David Sjostrand [david@ydsa.se] 

Datum: Sat, 5 Feb 2011 17:51+0100 

Dear Josef, 
dear Carl, 

I have worked a little bit more trying to find the locus that Josef is looking for. I think that it always is 
a line segment. I attach two files. 
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The first of David’s files was an appendix to his proof from page 5. 


Px=[aiA'\ 

Pi [^2 Ai\ 

— [( 2 ^ 5 ] 

Pn={\^yn) = {^-^)Pn-l+^Pn-i^n>2, 

We get two differenee equations 

x„ - (1 - m)x ^_2 - wx„_3 = 0 , Xj = , Xj = «2 Xj = (7) 

and 

- my „-2 - my„-3 , yi=h’ y 2 =bi and J3 = (8) 

(7) has the solution the general solution is x^ = Ar" + Br^" + C . Since < 1 and 1^2! < 1 we 
can conclude that limx„ = C . The initial condition Xj =a^, X 2 =a^ and X3 = gives us a 
possibility to find C. 

The expression of C is called x(al, a2, a3, m). You will find it in #10 in the attached dfw-file. 
x(bl, b2, b3, m) is the solution of (8) 

Therefore limP^ = [ x(al, a2, a3, m), x(bl, b2, b3, m)]. 

«->oo 


I feel very confident that converge if and only if 0 < m < 1 . 


Now we can plot for example [x(0, -5, 2, m), x(0, 3, 10, m)], 0 < m < 1. 
I have made this plot in Josef plot of 


VECTOR 


Here is the result 


ptsmZ 


0 ■ 


'S 


10 . 

. k, 50 

J 

J 

. k. 0, 1, 0.02 


JJ 

30 



Josefs red and my blue coincide quite well on a small interval. I think that the difference can 
be explained by thinking of that Josefs approximations for w-values close to 0 and 1 are too 
bad. 


The blue graph indicates that the graph is a line segment. This is true because if you simplify 


13 
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[x(0, -5, 2, m), x(0, 3, 10, m)], you receive , 
• . 13 

straight line y = - — x . 


3 


m + 2 


m + 


2 


which is the 


I think that [x(al, a2, a3, m), x(bl, b2, b3, m)] always is a straight line. 

Thus my belief is that the locus that Josef is looking for is a line segment for all m between 
0 and 1 . 

This is expression #10 and its simplification to give the straight line: 

#10: x(al, a2, a3, ni) 


3/2 

(v'(l - 4-in) - l)-(al-(l - 4-ni) •(2-iii - 1) - J(1 - 4-ni) • (al- (6-iii - 1) + 8-in-(a~ 


2 

32 -Ml •(:« + 2) 


2 + a3)) - 8-ni-(al-ni + a2 + a3)) 


#11: [x(0, -5, 2, 111 ), x(0, 3, 10, ni)] 


# 12 ; 


3 13 


III + 2 111 + 2 


Von: Josef Bohm [nojo.boehm@pgv.at] 

Gesendet: Sonntag, 6. Februar 201 1 1 1 :25 

Dear friends, 

This is an interesting dialogue, isn't it. 

Learning from David's approach using the difference equations I came to the same result 
today in the morning. 

My simplified function x(a1,a2,a3,m) is much more compact (seems to be identical, because 
the line segment is pretty the same). 

(expression #21 in the DERIVE file on page 13) 

I attach the respective file including my first thoughts (and failures!!) 

Best regards 

Josef 

(I skip printing some expressions of the next DERIVE session, this was presented earlier - 
on page 6.) 
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11 


#1: [Branch := Real, InputMode := Word] 

ptsni(1n1, ni, n, pt) := 

Prog 

pt := ini 
k :z 4 

#2 : Loop 

If k > n 
RETURN pt 

pt :z APPENDCpt, [n^ptiCk - 3) + (1 - m)^pti(k - 2)]) 
k :+ 1 

I take another set of initial points and want to approximate the locus of the "final points" 
(= limits of the sequence). My approximation is the set of the 50^^ points: 






0 

0 




VECTOR 


ptsni 


-5 

3 



, k, 0, 1, 0.02 





2 

10 

J 

J 



I added the slider (see also page 6) and proceeded: 



-4 -2 

My first impression: the locus is a closed curve! 

I add the set of the 100th points of the sequence on the same axes: 


ptsm 


0 0 
-5 3 

2 10 


, k, 100 


, k, 0, 1, 0.02 


100 



# 7 : 


VECTOR 
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This looks quite different; I get the impression of a very slow convergence. 
Compare the following results: 


Look at pts#30, 40 and 50: 


ptsm -5 3 , 0,02, 50 


2 10 


pts#98,99,100 


ptsm -5 


2 10 


[30, 40, 50] 


0 , 02 , 100 


,482470725 4,391027040 


i, 110177842 4,772135395 


807278865 5,082206674 


-1,976347000 5,932814058 


-1,003978754 6,928206876 


-1,956494305 5.953136842 


[98, 99, 100] 


And finally: 

pts#999, 1000 and 1001 



3 , 0,02, 1001 


-1,485148507 6,435643575 


-1,485148522 6,435643553 


-1,485148507 6,435643575 


[999, 1000, 1001] 


I try switching from numericai approximation to pure mathematics!! 

I follow David’s explanations about difference equations and try to generalize the problem. 
Given difference equation: x(n) = m-x(n-3) + (1-m)-x(n-2) 

x(n) -(1-m) x(n-2)-m x(n-3) = 0 = in general: x(n) + a x(n-l) + b x(n-2) + c x(n-3) 


General solution: x(n) = A r1 + B r2 +C r3 

with r1 , r2, r3 = roots of the characteristic equation: x^ + a x^ + b x + c. 


#11: SOLUTIONSCx - Cl - m)-x - m = 0, x) 


# 12 : 1 , - 


V(1 - 4.m) + 1 VCl - 4.m) - 1 


n f V'Cl - 4-m) + 1 y f JO- - 4-m) - 1 y 
#13: := a-1 + b- - + c- 


#14: m Real (0, 1) 


#15: lim K_Cn) = 1im 2 '(b-C- ^(1 - 4-m) - 1) + €■(-/(! - 4-m) - 1) ) + 


#16: m :e Real [0, 1/4] 

lim J<_Cn) = a 

# 17 : 
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I am going to find the coefficients now. 


cfs ;= CSOLUTIONS(x_(l) = al a = a2 a x_(3) = a3, [a, b, c])) 

#13: 1 

This simplifies to a huge and bulky expression with solutions for a, b, and c. To give an 
impression I print the first element of the solution vector which is a - and this should be the C 
in David’s paper. 

#19 ; ^ ;= 


5/2 3 2 2 3/2 4 3 - 

(1 - 4'm) ■(al*(2'm + 2'm - 1) - 2'a2'm '(111 + 2) + a3'(4'm - 2-m + 1)) + (1 - 4-m) ■(al-(8'm - 2-m + S-m- 


32 -m -(m +- 

2 2 3 2 2 

- S'tn + 1) + 2'a2'tn ■ (m + 2) • (1 - 4-m) + a3'(16-m - 12'm + 6-m - 1)) + 32'tn ■ (m + 2) ■ (1 - 4'm)'(al'm + a2 + a3- 


2 

2) -a - 4.m) 

) 

But to my big surprise when I selected the first element and simplified, then I received a very 
short and compact expression. 

cfs 

# 20 : 1 

al-m + a2 + aB 

# 21 : 

m + 2 

al-m + a2 + aB 

#22: co(al, a2, aB, m) := 

m + 2 

I defined a function depending on the first three coordinates and m and then tried this 
function on sequences from above. 


#23: [co(34. 

5, 100, 0.2), 

co(73, 3, 

-36, 0.2)] = 

■ 559 

92 ■ 


1 

1 

11 . 


#24: [co(34. 

5, 100, 0.2), 

co(73, 3, 

II 

1 1 

r\i 

0 

1 

[50,31318181, - 

3,363636363] 


Compare the next one with expression #10: 

#26: [co(0, -5, 2, 0.02), co(0, 3, 10, 0,02)] = [-1,435143514, 6,435643564] 

I add the line segment to the plot from above: 

#27: [co(0, -5, 2, t) , co(0, 3, 10, t)] 

This is the line in explicit form: 

X'(2'bl - b2 - b3) + al-(b2 + b3) - bl-(a2 + a3) 


#29 : y = 


2-al - a2 - a3 
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locusCirri, t) 


#B0: 


coCirri , itvi , ■irri , t) , coCirri , ‘irri , 1111 , t)' 

1,1 2,1 E,1 1,2 2,2 3,2 


#31: locus 



■ 0 

0 ■ 

> 


150 

650 " 


-5 


. 0.02 

_ 






101 

101 _ 




_ 2 

10 _ 

> 





#32: locus 


0 0 
-5 3 

2 10 


Function locus (ini ,m) returns the limit point to a given m. If you don’t enter a value for m 
then you can plot the locus in parameter form with 0 < t < 1 . 


VECTOR 


#33: 


VECTOR 


ptsm 


0 0 

-5 3 


LL 2 10 J JJ 


k. 3 


, k, 0, 1, 0,01 


j, 10, 100, 10 


The next graph shows #32 (the locus) and #33 which are the approximations of the locus 
starting with the locus of the 10*'^ points of the sequence (black) and the 100‘^ points of them 
varying m between 0 and 1 in 0.01 -steps (red): 
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Here are all approximations at together on the same axes (#32 and #33). 



Let’s have a jump into the 3'"'' dimension: 


VECTOR 


VECTOR 


ptsm 


0 0 0 
-5 3 4 

2 10 10 


k, 3 


k, 0, 1, 0,02 


I ^ \ ^ I 


j, 10, 100, 10 



0 


-5 
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Knowing the locus we can add it to our Tl-Nspire realization: 


m = 2 

I I 1^ 

0 . 1 . 


99. 57 


10 


34,73 


Mt)-- 


k1 ^+x2+x3 


v4^)-- 


t+2 

y1 z+y2+y3 


t+2 


(5.3) 


- 110.43 




10 


133 . 04 


( 50 * 809 ,- 3 - 36 ) 


( 100,-36) 


- 53 , 91 


The email communication went on: 

Von: David Sjostrand [david@ydsa.se] 

Gesendet: Sonntag, 6. Febmar 2011 23:21 

Hi, 

yes our correspondence is very interesting and stimulating. 

1. Josefs expression #21 proves that the locus always is a line segment. 

2. It is a little bit odd that Josefs expression #20 simplifies to the very compact and simple #21 when 
you don’t get this already in the first component in the vector in #19 and when I don’t get this simple 
expression at all. I have seen before that sometimes you have to let Derive simplify in two steps to get 
the simplest result. 

3. x(al,a2,a3) - #21 simplifies to 0 so they are identical. 

4. Josef s approximative curves are in a way more interesting than the correct line segments. 

5. It might be possible to go on with sequences of the type x(n)=a-x(n-l)H-b-x(n-2) + c-x(n-3), where 
a + b+ c=l. 

And here is David again: 

Von: David Sjostrand [david@ydsa.se] 

Gesendet: Sonntag, 13. Febmar 2011 21:04 

Dear Josef, 

I have spent some more time with convergent sequences defined by difference equations. 

Please, look at the attached files (convergent sequences.xls & convergent sequences.doc). 

Best regards, 

David 
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X_1 

7 

a_1 

0,14 

s_1 

0,14 


x_2 

2 

a_2 

0,12 

s_2 

0,26 


x_3 

-3 

a_3 

-0,23 

s_3 

0,03 


x_4 

1 

a_4 

0,16 

s_4 

0,19 


x_5 

5 

a_5 

0,81 

s_5 

1 



6,12 







5,4472 



S 

1,62 



3,941432 

3,396512 



b_1 

0,08642 

4,074074 


3,563348 



b_2 

0,160494 

4,074074 


4,033688 



b_3 

0,018519 

4,074074 


4,291805 



b_4 

0,117284 

4,074074 


4,261564 



b_5 

0,617284 

4,074074 


4,113921 





4,074074 


4,009922 

4,005838 





4,074074 


Von: Josef Bohm [nojo.boehm@pgv.at] 

Gesendet: Montag, 14. Februar2011 20:45 

Hi David, 

many thanks for this great idea. I am also working on some details and will come back to you 
as soon as possible. 

Best regards for the moment. 


Noor and Josef 









D 


-N-L#82 


David Sjostrand (and JB): The Proof 


p 19 


2-kI + + Z-kB) n f + 19-x2 - BO-kB) 

#11: ff(n, al, aZ, k1, kZ, kB) := ■! + |- 


17 


6S 


- 27-x2 + 14-:<B) 


74S 


S-JSS-L-Cl^-xl - 27-x2 + 14-kB) 


74S 


yCal +Z.al+4.aZ-B)+al-l 


7(al +Z-al+4.aZ-B)-al+l 


B.(ll'^l + 19-x2 - 30-x3) 


6S 


#1Z: ffCl, O.S, O.B, 1, B, -Z) = 1 


, O.B, O.B, 10 ul = 




1 B ■ 

#14: 

ptsn 

LTl 



.. -Z 6 . 

#1B: 

ffCn, 

O.B, O.B, 


7 B B7 ISB 17B7 BlBl BBB97 


1, 5, -Z, 


10 4 ZOO 400 4000 SOOO SOOOO 


- n/Z rr BB.]^1 9B.xZ 7B-:^B 


#1S: 


■K 


B4 B4 17 

-SIN n-ACOS + + + 

B74 J L L 4 JJJ 17 17 17 

#17: TABLE(ffC:<, 0,5, O.B, 1, 5, -2), , 10, 100, 10) 

10 0.41746 


cos n-ACOS 


[- 4 ]] 


B.^BB.(Z7.:<Z - 14.:^B) 


B74 


-BS 

20 0.41176 + 2.4676.10 ■ i. 


-42 


BO 0.41176 + 4.SS2-10 

40 0.41176 

50 0.41176 - 5.494S-10 

60 0.41176 - 4.9B6-10 

70 0.41176 - 5.4BS-10 


-49 


-52 


-55 


-59 

SO 0.41176 - 7.1662-10 ■ i. 

-62 

90 0.41176 + 6. SB56-10 ■ i. 

-66 

100 0.41176 + B. 4224-10 ■ i. 


#15: 


2-xl + B- (x2 + 2-x3) 


17 


2-1 + B- CB + 2- C-2)) 


#20: 

7 

# 21 : 

17 

#22: 0.41176 


17 



Rule from page 17 
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Von: Yvonne Fogelander- Sjostrand [yvonne@ydsa.se] 

Gesendet: Mittwoch, 2. Marz 2011 21:46 

Dear Josef, 

thank you very much. I will look at your files with great interest. 

You are right. It should be n>=0. 

Best regards, 

David 


Von: David Sjostrand [mailto:david@ydsa.se] 

Gesendet: Mittwoch, 9. Marz 2011 21:49 

Dear Josef, 

now I have been investigating Davids Sequence S.dfw and I think I understand what is happening. 

In Davids Sequence_4.dfw I have investigated some more expressions. 

#30 is a little bit more symbolic and its simplification confirms my rule a little bit more. 

The function ptsn(i ni , al , a2 , n , pt) is very nice it makes it easy to investigate sequences of 
order 3. Would it be possible to make it more general so that you can use it to investigate sequences of 
arbitrary order? 

I will soon spend some time with Davids Sequence_2.dfw 
Best regards, 

David 


#30: 




ptsn 





xl 

x2 

x3 




0.5, 0.3, 50 


[4S, 49, 50] 


2 5 10 

#31: xl + x2 + x3 

17 17 17 


#32: 


0. 11764 -xl 
0. 11763 -xl 
0. 11763 -xl 


0. 29410- x2 + 0.5SS21-X3 
0.29409-X2 + 0.5SS1S-X3 


0.29409-X2 + 0.5SS1S-X3 J 


0.11764-xl + 0.29411-X2 + 0.5SS23-x3 


Von: Josef Bohm [nojo.boehm@pgv.at] 

Gesendet: Donnerstag, 17. Marz 2011 20:58 

Dear David, 

Many thanks for your new dfw-file. 

I will try to "compose" a nice article on this topic for the next DNL. 

Best regards 
Josef 

I leave the challenge to prove - or just double check - David’s conjecture from 
page 17 for you! (Or you go to page ...) 
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Ecoie de technol(^ie superieure 

Milt 145-03-115-07 : Culcul differentiel et integral Session A-2010 

Michel lleuiidin 


Devoir #2 

Corrige siir 300 ]>oiiits, ]K>iide ration : 15% 

Dofuie (pour tes 3 groupes) te imdi 8 novembre 2010. 

A remettre an SEG (8-2500) avant J7h, (e vendredi 3 decemhre20l0. 

A faire en iH|iii]H' de 2 a 4 iHiidUmLs. 

Note 1 : cles solutions de piusieurs exercices de ce devoir seront donnees en T.P. les jeudi 9 ^AM. 
gr. 03; PM. gr. 07) et vendredi JO decemhre )AM, gr. 05). Les devoirs corriges seront egalemeru 
remis aux etudiaiits pendant ces T.P. Notez, de plus, que la matieie vue aii dernier cours 
(semaine du 6 decembrej N’EST PAS evaluee dans ce devoir mais fait partie de rexanien final. 

Note 2; vous n’avez pas a calciiler a la main les integrales des problemes suivants; #.3, #4, #12, 
#13 et #14. Dans ce devoir, utilisez « Display Digits: Hoat 6» et «« Complex Format; 
Rectangular ». Due (autiej table d'integrales est annexe a la fin de ce document. ; avec plus de 
symetrie entre les differentes formules, qui contient des primitives pour les racines carrees 
d' expressions quadratiques ati nuirterateur et qui esten consequence directe avec le problenie 3B) 
du devoir I. Je vous Toffre comme cadeau de Noel (nierci a mon ami Albert Richet a RUBE). 

l^rohlenie 1 |15 points! Sotntiu' de gniiche, soiiime de droite Considerons Pintegrale 
s I 

definie suivante; / = |v 1 + f ^ dx . 

I 

11 ) [)onnez une raison pour laquelle Pexpression a integrerest croissante sur Pintervalle [!, S]. 

b) Quelle valeur, pour /, est retourn^ par votre calciilatrice? 

c) Remptissez le tableau suivant (utilisez vos fonctions « gaucEve », « droite » ou encore 
Papplication Flash Cafcufus Tools) pour les calciils. Notez que vous avez bien les iiiegalites 
Gauche(n) <I < Droite(/i). pout chacune des valeurs de n dans le tableau. 


n 

Gauche(n) 

DroiLLifn) 

5 



20 



.35 




Problenie 2 |20 points} Construction de primitive (suite du probleme 1). En mode exact, 
votre calculatrice ne pent trouver la valeur de / (certains logiciels le font mais la reponse est fort 
compliqu^ et incomprehensible pour un etudiant du premier cycle). Posons 

toto(.r)= [ , dt . 

0 Vl + f^ 

Definissez cette fonction sur votie Tl. Vous allez trouver la valeur de / en terntes de « toto ». 

11 ) En repondant en fonction de « toto », quelle est la valeur de f , dt ? 


1 
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b) Demontrer la formule suivante, simptement en clerivant k cote droit : 




l(fl 


c) Motitrez main ten ant que I - 


I 

np + 1 up + 1 

48>^-2V2 2 




5 


+ ^{toto(S)- toto(])). Ainsi, si la fonction ‘1oto’’ 


existait dans votre calculatiioe, Tintegrale/ aurail ete calculee de fa^on exacte 


Prohleine 3 [ 15 points | Applications plij siqnes de rintegrale definie Un ohjet se de place le 
long de l"axe des .? et sa vitesse, a Tinstant t, est dcmnee par v(/) = e~‘ +fcosr + 2sinf (en ni^s). 
Au debut des observations, a f = 0. la particule est a la position x - 3 m. Comhien de fois. au 
cours des 20 premieres secondes, la particule traversera-t-elle Torigine? 

Probleme 4 |15 points! Application diverse de I’ integrale definie Le nombie d'heures //oil 
il fait jour a Madrid en fonction de la date est calcule approximativement par la formule 


// (r J = 1 2 + 2.4 si ntO.O ! 72f# - 80)) 


ou t est le noinbre de jours depuis le debut de Tannee. Troiivez le nombie d'heures moyen oit il 
fait jour a Madrid ; 

a) en janvier ; b} en juin ; e) durant une ann^ complete. 

Probleme 5 [15 points] Tbeoreme fonda mental dii calcnl integral Soit la fonction 
derivee f\x) donnw par le graphe suivant ; 



Si ytO) = 2, remplissez le tableau suivant en utilisant le theoreme fondamental du calcul integral et 
des calculs d’aiie seulement (vous NE devez pas trouver/tjc) explicitement!). 


X 

0 

I 

2 

3 

4 

5 

6 

fixi 

2 
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Probleiiie 6 |25 pointsl Cnlciil d’mteKi'oles du tvpe [ — . Soient A, Bet C trois 

■*A.r + BA+C 

nomhres r^ls avec A non nul. Pour calculer f — , on complete le cane et utilise un 

A?- + Bf + C 

chaneeinetu de variable. Cela nous ameive a I'un des 3 cas conniis suivanis : 


rdy ] _ 

J V' ¥ 


^ = — arctan — (ti>0); ^ = — arctann — . (n > 0). 

£i‘ + ¥‘ ^ a ) ^ - v‘ a \ a ) 


Appliquei cela pour calculer les iiitegrales suivantes : 
r dx , r dx 


ax 

-} 

bl J 

X + 20a' + 1 09 


dx 

(i) r 

-.X- + 6.X-5 

J 


J- + 10I + 21 

dx 

2x--x + ti 


r cir 
■' x‘ + ]0x + 25 


Probleiiie 7 |35 points] CaJcxil d^intejirales dii type f — cas 

A-X ~ + Bx + C 

d'intecrales du tvpe f — tiotre premier reflexe devi ait etre de verifier si la dei ivee 
' ■'Ai- + Ba+C 

du (Jenominateur n’est pas, a facteur pres, au aumerateur, ce qui donnei a, comme reponse, un 
seui logaiithnie. Si tel n’est pas le cas, on regarde les l acines du denominateiir. On fera des 
fractions partielles si les racines sont reelles (« expand » peut etre utilise mais ici c’est tres rapide 
a la main) et on completera le cane si elles sont complexes. Appliquei ces piincipes pour le 
calcu! des integrates suivantes : 

r {.x + 10)i7.x , r lit r (.3-2.t)iii 

a) ^ b> ^ Cl 

. 1 - +20.X+ 109 .r + 10.X + 2 [ .X- + I0.X + 25 


.1+ 20.x +109 
r 5.x dx 

^ -.X-+6.X-5 


.X- + 10.X + 21 
r (7.x-]jrf.x 

2.x --.x + 6 


Problenie 8 [50 points] Integration par diangeinent de variable Chacmie des integi ales 
suivantes peut se faire par un changement de variables. Calculei-les en utilisant votre table 
d'integrales [ou celte a la derniere page de ce document). 

I \/2 , " Jr)6 


l2x^]- x- dx b) J 


rdx c} 


0 VI -A- 


:} J s i n 1 3.x ) CO s( 3.x ) dx d } J .x ‘ [ 4 - 5.x ■ J dx 


f ^\2i) 

e) j ti j-JL^ R} j(Wx-+lx-'^)'‘dx 

i(^l + .x-j 

Probleine 9 ] 15 points] Integration par parties Calculei les integrates suivantes ; 


b} ( 1 0.X" + 1 ) cos[2.x) dx 


.i 
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Problenu' 10 [10 i^oiiit^l Integration par fractions partielles Calculez les integrates 
suiv antes, apres avoir uliiise la commaiide « expand » de votre calcuiaLrice : 


r IV. 1 

- dx 

+ 

1 

i-j 

1 

1 ^ 


'x-+9] 


Prohli'ine 11 [15 points} Integration tPe\pressions (jna(lratir]nes sons un radical Calculez, 
apres avoir complele le carre, les integrates suiv antes : 



Problenie 12 [15 |)oints1 Integrales impropres tndiquez pourquoi tes integrates suivatites 
sont impropres, dites si ettes convergent et si oui, trouvez tear vateiir. 




dx 


c) 

0 (COSA-)^- 


Problenie 13 [15 points} Applications geonietriqnes de I’integrale definie Dessinez ta 
region detimit^ par les coiirbcs = et y = Siru.r + 1) (voiis pouvez reproduine te 

graphlqne). 

a) Trouvez I’aire de cette region. 

b) Trouvez la longueur de la frontiere de cede region. 

0 } Trouvez le volume dii solide de revolution engendre par la rotation de la region autour de 
Taxe des .r. 


Problenie 14 }1? points) Applications geonietriqnes de I’integrale definie Dessinez la 
region delimitee par les courbes y - 5 siii>/7 et y - . 

a) Trouvez I’aire de cette region. 

b} Trouvez le voiume tlu soiide de revolution engeiidre par la rotation de la legion autour de 
I’ axe des x. 

c) Trouvez le volume du soiide fie revolution engendre par la rotation de la region autour de la 
droite horizontale y= I. 


ProbleniL' 15 }25 points} Series de puissances, SH'ies de Taylor 

A) Trouvez rintervalle de convergence de chacuiie des series de puissances suiv antes : 

^ H - 1 n" 

^ n—\ “ n-\ ^ ^ 

P) Trouvez les developpements de MacLaurin des fonctions suivantes jusqu’a I’ordi^ 4 et 
indiquez aussi rintei'vulle de convergence. Allez-y tie 2 famous differentes: en utilisunt la 


definition de seiie tie MacLaurin 



eten ulilisant les series tie base de la table. 


1 

u ; 

lO+.r 


ii} — 

9-.t‘ 


iii} cos 



4 
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Uiie (autre) table d^integrales 

Note ; In constante d'integration aete omiseet a designe an nombre reel positif ta> 0) 

« Les formubs de base » 

r 1 

[. \x‘’dx = - — -^/A- = ltiA- 


4. [sin A' £/a‘ = -cos.y 


. r I J I f 

6. ^ ax = — arctan — 

‘'i7‘+A“ 


arc tan h 


Le lestecles tonctions trii:Linometrit]UL:s de base 


8. J lan A dx = - In [cos a) 

9. J cot A dx = ln(sio .a) 

[0. 1 sec A tir — ai'ct an h t si n a ) 

1 ] . 1 CSC A dx - - ai ctanlHcos Aj 

12. J sec* A dx= ran A 

13. jese* A cif = -COfA 


se c A t aji A dx = se c a i 5 . esc a co t a dx = - esc a 


Les sommes et differences de carres sous un radical » 


.JA = arcsin' — 


= arcsinh 


f£l" + A 


[ 




\ f = dx = ardiinh 

-d-a- +A" 


f-a~ + A" 
A 


[9- - A" Ja = — +— f— i= 

J 0 0 ^ rr 


20 . l^la-+x- dx = ^^|ir+x- +^J I J , 
^ ^ \a- +x‘ 


£J" +A‘ dx 





















p 26 


Dietmar Oertel: Structured Combinatorics 


D-N-L#82 


Strukturierte Kombinatorik - Structured Combinatorics 

Laws and Structures in Combinatorics 
and 

“Secrets of the Prime numbers” 

Dietmar Oertel, Miinchen, Germany 



Abstract 

The following paper deals with laws, structures and patterns in combinatorics in 
co nn ection with the Triangle of Pascal. In order to work meaningful with 
combinations of elements they must be given a structure. Having introduced a 
structure, special combinations can be addressed and if necessary supplied with 
attributes and then worked up under particular points of view. 


Important advantages of a structure are: 

• a clear order of an arbitrary assembled group of elements, 

• a unique address of combinations, 

• the order of structured combinatorics gives an answer to the old question 
of a pattern of the prime numbers. 

The paper consists of two parts. 

The first part, Structured Combinatorics,, demonstrates the patterns appearing in the 
Triangle of Pascal. 

The second part, Structured Prime Numbers, deals with the ‘phenomenon’ of the prime 
numbers, which can be described based on the Structured Combinatorics. 
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1 Structured Combinatorics 

Unsorted and unstructured elements of a set cannot used for follow-up calculations, 
observations, Investigations, and simulations in most cases. This makes necessary to find a 
common structure for all possible combinations of elements. 

The number of all possible combinations is given by 

n ! 

KmbCn, k) 

k ! ■ (n - k) ! 

which is in DERIVE available as COMB(n , k) . 

The Triangle of Pascal is easily generated with 
PascalTriangleCn_from , n_to) := VECTORC[VECTORCKmb(j , k) , k, 0, j)], j, n_from, n_to) 
For example: 

■ 0 [ 1 ] 

1 [ 1 , 1 ] 

2 [ 1 , 2 , 1 ] 

3 [1, 3, 3, 1] 

PascalTriangleCO , 7) = 

4 [1, 4, 6, 4, 1] 

5 [1, 5, 10, 10, 5, 1] 

6 [1, 6, 15, 20, 15, 6, 1] 

.7 [1, 7, 21, 35, 35, 21, 7, 1] . 

Row 5 gives the number of groups of selections of a given set of 4 distinct elements (without 
consideration of the order): There is one selection with no elements, there are 4 ways to select 
one element, 6 ways to select 2 elements, again 4 ways to select 3 elements and only one way 
to select all four elements. 

We can see the groups of all different numbers of elements using the next function (base set 
of 5 elements): 


KmbObjAl lNor([a, , 5, 10, b]) 
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We can also find all selections with a given number of elements (e.g. for 4 elements): 

ax 5 10 


KmbObjGrpC[a, x, 5, 10, b], 4) = 


ax S 
a X 10 
a 5 10 

X 5 10 


Earlier we spoke about finding a unique address for eaeh single eombination. 

Funetion KmbPosOf ([Combi nati on] , [base set]) returns the “address” of the eombination 
- whieh is the position in the Triangle of Pascal. 

Funetion KmbNteObjIn(pos,grp, [base set]) finds the eombination at the given position. 
Let’s give two examples (applying DERIVE’s DISPLAY-eommand): 


What is the position of [10, b, a] in all eombinations of [a,x,5,10,b]? 

#1S: KmbPos0fKmbC[10, b, a], [a, x, S, 10, b]) 

requested (sorted): 

[10, a, b] 

in base set (sorted): 

[S, 10, a, b, x] 

Group: 3 
Group size: 10 
Position : 


#19: 7 

The function gives an answer: It is in the 7* row of the column of three elements’ 
combinations. Is this true? 


We double check by generating all groups of possible combinations of the given 5-element 
set: 


# 20 : 


KmbObjAl lNor([5 , 10, a, b, x]) 



5 10 a b 

5 10 a X 

5 10 b X 

5 a b X 
10 a b X 


[[5, 10, a, b, x]] 


# 21 : 



D 
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Yes, indeed. The eombination of 10, a, and b is placed in the seventh row of the group built 
by the 3-element combination. But we could also use the inverse function (without regarding 
the order of the elements!) 

#22: KmbNteKmbIn(7, 3, [5, x, a, b, 10]) 

Nth Element : 

7 

in Group: 3 

in base set (sorted): 

[5, 10, a, b, x] 

Group Size: 10 
Combination: 


[1Q| a, b] 


Another example with more elements: 

#24: KmbPosOfKmbCLlS. 11, 1, 2], [1 20]) 

requested (sorted): 

[1. 2. 11. 15] 

"in base set (sorted): 

[1. 2, 4, 5, 6, 7, S, 9, 10, 11, 12, 1^, 14, 15, 16, 17, IS, 19, 20] 

Group: 4 

Group size: 4S45 

Position : 

#25: HE 

Would you like to check this manually? 1 prefer asking DERIVE doing the job: 

#26: KmbNteKmbIn(112, 4, [1 20]) 

Nth Element : 

112 

in Group: 4 

in base set (sorted): 

[1. 2, S, 4, 5, 6, 7, S, 9, 10, 11, 12, IB, 14, 15, 16, 17, IS, 19, 20] 
Group Size: 4S45 


Combination : 
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KmbPos0fKmb([2,3,4],[l 10]) = 37 

KmbNteKmbIn(37,3,[l, . . . ,10]) = [2,3,4] 

KmbPosOfKmb([ll,15,12,17],[l 40]) = 64040 (11 sec) 

KmbNteKmbIn(64040,4,[l, . . . ,40]) = [11,12,15,17] 

KmbPos0fKmb([l,2,7,16,20,33], [1 33]) = 15757 (154 sec) 

KmbNteKmbIn(15757,6,[l 33]) = [1,2,7,16,20,33] 

I am no mathematician. I did not come across a way to describe this structure in a short way 
until now. So I am using for my own the following notation: 


Sub = 4 

Posx = Pj^ 33J ([11, 12, 15, 27]) = 3211 6 

Main=33 
Sub = 4 

Kmbx = 33 J (321 1 6) = [1 1 , 1 2, 1 5, 27] 

Main = 33 

Remark: If there is any official notation then please let me know. I would be very grateful 
when a professional mathematician would indicate a more appropriate notation, which could 
be applied for mathematical problems. 


TIME 2010 - The Proceedings 


Dear colleagues, 

Now you can download the CD version of Proceedings of TIME 2010 and the covers for the 
CD box. In order to do it, go to TIME 2010 website (http://www.time20 1 0.uma.es/) , select the 
Proceedings section and then visit the section CD version of Proceedings. 

You will find two different options to download: 

• PROCTIME2010.iso: This is an .iso-version which can be burnt with any software 
which supports iso files. 

• PROCTIME2010.zip: A zip file with the contents of the Proceedings. 

If the CD does not auto-run, open the file index.html to see the contents of proceedings. 


Best regards and we hope to meet you once again in Estonia for TIME 2012. 
Jose Luis, Gabriel and Pedro 
Editors of TIME 2010 Proceedings 
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# 1 : 


ptsn(ini, coeffs, n, pt, pt_, d, k, k_) := 

Prog 

pt := -ini 
d := DIHCini) 
k := d + 1 

coeffs := APPENDCcoeffs , [1 - I(coeffs)]) 

Loop 

If k > n 
RETURN pt 
k_ :z DIM(pt) 

pt_ := ICcoeff^i 1 ■ ptj, (k_ - d + i), i, 1, d) 
pt := APPENDCpt, [pt_]) 
k :+ 1 


# 2 : 


r 

f- 

ptsn 





#3: 


ptsn 


M 7^ 
100 -56 


7 

2 

-5 

1 

5 


[0.2, O.S], 100 


50 . 51515151 -S. 565656565 

50 . 51515151 -S. 565656565 
50 . SlSlSlSl -S. 565656565 


[0.14, 0.12, -0.25, 0.16], 7 


[9S, 99, 100] 

7 
2 
-5 
1 
5 

6.12 
5.4472 


#4: 


7 
2 

ptsn -5 
1 
5 

David's formula for the limit 


[0.14, 0.12, -0.25, 0.16], 100 


[4.074074074] 


100 


#5: 


seq1im(ini, coeffs, d, den) := 

Prog 

d := DIHCim) 

coeffs := APPENDCcoeffs , [1 - I(coeffs)]) 
den := laCcoeffs^k, k, 1, 1), 1, 1, d) 
ICimi 1 -ICcoeffsi. 1 , 1, 1, i)/den, i, 1, d) 


#6: seqlim 


#7: seqlim 


54 75 

5 5 

100 -56 

7 


[0.2, O.S] 


559 


L 11 


92 


11 J 


2 
-5 
1 
5 J 


[0.14, 0.12, -0.25, 0.16] 


= [4.074074074] 
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A Mathematical Model for Snail Shells (2) 

Piotr Trebisz, Germany 



Im 1. Teil meiner Artikelreihe babe ich eine Methode vorgestellt wie man Schnecken- 
hauser als parametrisierbare Oberflache eines Torus modellieren kann, der auf einer 
planaren logarithmischen Spiral-Kurve basiert. 

Diese sind aber nur ein Spezialfall. Deshalb werde ich im 2. Teil meiner Artikelreihe 
zeigen wie man den interessanteren, aber auch komplexeren, Typ modellieren kann, der 
auf konischen logarithmischen Spiral-Raumkurven basiert. 

Ich nenne dieses konische Schneckenhaus "Trebisz-Spirale - Typ 2". Da die Konstruk- 
tion von Typ 2 analog ist zur Konstruktion von Typ 1, werden wir einige der Dinge, die 
wir im 1. Teil gelemt haben, auch hier anwenden konnen. 

In the first part of my series I presented a method how to model snail shells as 
parameterized surface of a torus based on a logarithmic spiral. 

As these surfaces are a special case I’d like to demonstrate how to model a 
more interesting but more complex type which is based on conical logarithmic 
spirals in space. 

I will call this conical snail shell “Trebisz Spiral Type 2”. Construction of type 2 
surfaces is analogous to construction of type 1 bodies, so we can apply some 
facts which we have learnt in part 1 . 
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Als Geriist fiir dieses Schneckenhaus dient eine konische logarithmische Spiral-Raumkurve die nach 
ihrer Lange parametrisiert ist und die sich um die X-Achse windet. Der Parameter t steht hierbei fiir 
die Lange der Spirale entlang der X-Achse, b ist die Basis der Logarithmus und m ist das Verhaltnis 
zwischen der Lange des Konus und seinem Radius. 


As base for this type of snail shells serves a conical logarithmic 3D-spiral which is parameter- 
ized by its length winding around the x-axis. Parameter t describes the length of the spiral 
along the x-axis. b is the base of the logarithm and m is the ratio between the length of the 
cone and its radius. 


#9: SNAILCb, m, t) [t, m.t.COSCZ.jr.LOGCt, b)), m.t.SINCZ.n.LOGCt, b))] 

Example: SNAlL(e, 1/2, t) with0<f<1 



Dass diese Spirale parametrisiert ist nach der Lange entlang der X-Achse ist offensichtlich. Dass die 
Spirale auch nach ihrer Bogenlange parametrisiert ist, lasst sich dadurch zeigen, dass die Tangente an 
jedem Punkt der Spiral-Raumkurve eine konstante Lange hat. 

It is obvious that this spiral is parameterized by its length along the x-axis. We can show that 
it is also parameterized by its arc length: the tangent in every point of the spiral has constant 
length. 

2 2 2 2 
+ l)'LN(b) + 4'7 t '01 ) 

LNCb) 


# 10 : 


dt 


SNAILCb, m, t) 


Wie man sieht ist die Spirale auch nach ihrer Bogenlange parametrisiert, aber skaliert um einen kon- 
stanten MaBstabs-Faktor. 

Um den Mantel des Schneckenhauses zu konstruieren benotigen wir das begleitende Dreibein, wel- 
ches wir mittels der Frenetschen Formeln ermitteln konnen. 


As you can see the spiral is parameterized by its arc length but scaled by a constant 
measurement factor. For constructing the surface of the snail shell we are needing the 
accompanying three pod which can be obtained by Frenet’s formulae. 


#11: X_AchseCb, m, t) SIGN 

Td '] 

— SNAILCb, m, t) 

X 

^d '' 

2 ^ 
SNAILCb, m, t) 


.^dt j 


^dt> 



#12: Y^chseCb, m, O ;= - SIGN 


ffd 'l 

2 

— 

SNAILCb, m, t) 

^Ldtj 



#13: 


Z^AchseCb, m, 


rd 


t) ;= SIGN 


— SNAILCb, m, t) 
^dt 


y 
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Dieses lokale mitdrehende Koordinatensystem hat dieselbe Handigkeit wie das globale Koordinaten- 
system, in dem die Spirale definiert ist. 

This local co-rotating system of coordinates has the same direction as the global system 
where the spiral is defined. 

#14: X_Achse(b, m, t) x Y_Achse(b, m, t) = Z_Achse(b, m, t) 

Analog zum Typ 1 soil die Mantelflache einen kreisformigen Querschnitt haben. Zusatzlich soil dieser 
Kreis an jedem Punkt der Spirale die X-Achse beriihren damit um die X-Achse herum keine Liicke 
entsteht. Es stellt sich also die Frage wie der Radius des Kreises von t abhangen soil, damit diese Be- 
dingung erfullt ist. 

Der Ansatz ist der, dass wir uns zu jedem Punkt auf der Spirale einen korrespondierenden Punkt auf 
der X-Achse suchen, so dass die Verbindungslinie zwischen den beiden Punkten senkrecht zur Tan- 
gente steht. 

In analogy to type 1 the surface shall have a circle as cross section. Additionally this circle 
shall touch the x-axis in each point of the spiral such that there will be no gap around the x- 
axis. How shall the radius of the circle depend on f to fulfill this condition? 

We try to find for every spiral point a corresponding point on the x-axis such that the 
connecting line between these points is perpendicular to the tangent. 

#15: LINIECb, m. t, a) := [a, 0. 0] - SNAILCb, m, t) 


Gesucht ist ein a fiir das gilt: We are searching for an a with: 

d 

#16: LINIECb, m, t, a) • — SNAIL(b. m, t) = 0 

dt 

Wir losen die Gleichung nach a / We solve for a 
2 

#17 : a = m ■ t + t 

Wir setzen a in LINIE(b , m , t , a) ein und bestimmen die Lange der Verbindungslinie. Damit erhalten 
wir den passenden Radius fur den Mantelkreis. 

We substitute for a in LlNlE(b , m , t , a) and calculate the length of the connecting line 
segment. This gives the appropriate radius for the circle of the surface. 


#1S: 


2 

LINIECb, m, t, m -t + t) 


2 

+ 1 ) 


Wie man siebt wacbst der Radius des Mantels, analog zum Typ 1, linear mit der Lange der Spirale. 

We see that the radius of the surface - like type 1 - is increasing linearly with the length of 
the spiral. 

2 

#19; MANTEL (b, m, t, tp) := + 1) ■ t- (COSCtp) ■XJchseCb , m, t) + SINCtp) ■ Y_AchseCb , m, t)) 

Damit baben wir alle Zutaten fiir unser Scbneckenbaus. Es ist die Summe aus Mantel und Spirale. 

So we have now all ingredients for our snail shell. It is the sum of surface and spiral. 

#20: Spi ralTorusX_LCb , m, t, cp) := SNAIL(b, m, t) + MANTEL(b, m, t, cp) 



Wenn wir das Schneckenhaus stattdessen nach der Anzahl der Umdrehungen parametrisieren wollen, 
rniissen wir t durch ersetzen. 

If we want to parameterize by the number of rotations we have to replace t by b\ 

t 

#22: Spi ralTorusX_UCb , m, t, (p) 5piralTorusX_L(b, m, b , tp) 

Wie zu erwarten war gibt es auch beim konischen Schneckenhaus das selbe Problem mit dem Parame- 
ter m. 1st m zu groB, dann durchdringt sich die Mantelflache bei jeder Umdrehung selbst. 

How expected we are facing now again the problem with the parameter m. If m is too big 
then the surface will interpenetrate. 

Example: 

Spi ra1TorusX_L(e, 2/3, t, (p) 

with 

0< t<l/\-n<(p<n 



1st m hingegen zu klein, dann klafft eine hassliche Liicke zwischen den Windungen des Spiraltorus. 

If m is too small, then we have an ugly gap between the windings of the surface. 


Example: 

Spi ralTorusX_L(e, 1/4, t, 0) 
with 

0< t<lA-TT<(p<TT 



Wir wiirden geme m so bestimmen, dass die Mantelflache sich bei jeder Umdrehung selbst beriihrt 
aber ohne sich selbst zu durchdringen. 

We would like to determine m such that the surface is osculating itself at every winding 
without interpenetrating itself. 


Example: 

Spi ralTorusX_L(e, 571877/1286241, t,c|)) 
with 

0< t<lA-TT<(p<TT 
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Es iiberrascht wohl kaum, dass der Ansatz zur Losung dieses Problems analog zur Losung des Selbst- 
durchdringungsproblems beim Schneckenhaus von Typ 1 ist; Wir bestimmen den Schnitt des Schne- 
ckenhauses mit der XY-Ebene. Die Z-Koordinate muss also 0 sein. Doch anders als bei Typ 1 schnei- 
det Typ 2 die XY-Ebene immer und nicht nur bei jeder halben Umdrehung. Deshalb muss man die Z- 
Koordinate bier nicht nach t sondem nach (p auflosen. 

Im Prinzip konnte man dabei genauso vorgehen wie beim Schneckenhaus vom Typ 1 und die Z-Koor- 
dinate nach (p auflosen, allerdings ist diese beim Typ 2 ein ziemlich komplizierter Term und die Lo- 
sung die wir dabei erhalten ist sehr unhandlich. Um uns die Arbeit zu erleichtem suchen wir uns statt- 
dessen eine alternative, aber aquivalente, Parametrisierung fur das Schneckenhaus, dessen Z-Koordi- 
nate sich einfacher nach cp auflosen lasst. 

Wir haben ja vorhin die Lange der Verbindungslinie zwischen der Spirale und der X-Achse bestimmt 
um den Radius des Kreises zu berechnen. Genau diese Verbindungslinie konnen wir stattdessen als 
lokale mitdrehende X-Achse verwenden, denn sie hat ja bereits die passende Lange und steht senk- 
recht zur Tangente. 

It is not surprising that the solution of this problem is very similar to the solution of the 
interpenetrating problem of the type 1 snail shells. We determine the intersection of the shell 
with the xy-plane - hence the z-coordinate must be 0. But in a different way from type 1 
type 2 intersects the xy-plane always and not only at every half rotation. As consequence we 
have to solve the z-coordinate not for t but for cp. 

Basically one could proceed like with type 1 shells and solve the z-coordinate for cp but this 
expression is here a pretty complex term and the solution a bulky expression. In order to 
make work easier we are looking for an alternative but equivalent parameterization which 
enables a more comfortable solution for cp. 

We have calculated the length of the segment between the spiral and the x-axis in order to 
calculate the radius of the circle. This line can be used as local co-rotating x-axis because it 
has the appropriate length and is perpendicular to the tangent. 

2 

#Z4: X_Achse_NeLj(b , m, t) := LINIECb, m, t, m -t + t) 

#25: y_AchsG_NGLiCb , m, t) := IlINIECb, m, t, m -t + t)|-SIGN|^|^ — SNAIL(b, m, t)j x X_AchsG_NGLiCb , m, t)j 

#26: MANTEL_NEU(b, m, t, ct:) := CO5(0) ■X_AcbsG_NGuCb , m, t) + SINCcp) ■ Y_AchsG_NGj(b , m, t) 

#27: ScbnGckG_L(b, m, t, ip) := SNAILCb, m, t) + MANTEL_NEU(b , m, t, p) 


Um das Rechnen zu erleichtem, arbeiten wir lieber mit dem Schneckenhaus, das nach der Anzahl der 
Umdrehungen parameterisiert ist. 

In order to make calculation easier we prefer working with the snail shell parameterized by 
the number of rotations. 

t 

#28; Scfinecke_UCb , m, t, 0) := Schnecke_L(b , m, b , cfi) 

Diese Funktion beschreibt das selbe Schneckenhaus wie die urspriingliche Funktion, die Terme sind 
aber einfacher und handlicher. Die beiden Funktionen unterscheiden sich nur in einer konstanten 
Differenz des Winkels (p. 

Auch die Z-Koordinate lasst sich jetzt einfacher nach (p auflosen. Es empfiehlt sich aber dies per Hand 
durchzufuhren, da die Ergebnisse, die DERIVE liefert, unnotig kompliziert sind (#31). 
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Both functions are describing the same surface. They only differ in a constant difference of 
the angle cp. 


The z-coordinate can be solved easier for cp, too. I recommend performing this manually 
because the DERIVE results are unnecessarily complicated (#31). 


#30: 0 




t 

m- b ■ 1 


2 

Z-JT-m -SINCtts) 


2 2 2 2 
yC(m + l)'LNCb) + 4'JT 'm ) 


- COSCtts) 


■ SINCZ.JT.t) 


2 

Cm + l).LNCb).SINCtt3)-COSC2.JT.t) 


2 2 2 2 
7C(m + l)'LNCb) + 4'IT -m ) 


" 2 2 " 
(m + 1) ■ LN(b) ■ C0T(2 'JT' t) - Z-jr-m 


#31: 


$ = 2-ATAN 


2 


2 2 


L + l)'LNCb) + 4-jt -m ) j 

Wenn wir jetzt diese Losung einsetzen verschwindet die Z-Koordinate und wir erhalten die 
Schnittkurve des Schneckenhauses mit der XY-Ebene. 


Substituting this solution the z-coordinate disappears and we obtain the intersection curve of 
the snail shell with the xy-plane. 


#32: SCHNITTCb, m, t) := 


2.LN(b).m.b 


2 t 

m .b . 


2.SIN(2.7r.t).LN(b).(LN(b).SIN(2.7r.t) + 2.7r.COS(2-7r. t)) 


2 2 2 2 
L (m .COS(2.7T.t) + l).LN(b) + 4.7r.m .5IN(2.7T.t). (7r.SIN(2.7r.t) - LN(b).C0S(2 -n. t)) 

2 2 
(m + l).LN(b).COS(2.7T.t) - 2.7r.m .SIN(2.7T.t) 


2 2 2 2 

(m .COS(2.7T.t) + l).LN(b) + 4.7r.m .SIN(2.7T.t). (7T.SIN(2.7r.t) - LN(b).COS(2.7T.t)) 


Man erhalt exakt die selbe Schnittkurve wenn man die Z-Koordinate der originalen Parametrisierung 
Spi ralTorusX_U(b,m, t ,(()) nach (p auflost und die Losung dann fur (p einsetzt. Ich werde das hier 
aber nicht extra beweisen, stattdessen ermutige ich die interessierten Leser sich selbst daran zu ver- 
suchen. 

Die Wahl dieser alternativen Darstellung fiir das Schneckenhaus hat, neben einfacheren Termen, noch 
einen anderen groBen Vorteil. Und zwar lassen sich damit geeignete Startwerte fiir das Newton- 
Verfahren bestimmen, auf das ich spater noch im Verlauf meines Artikels naher eingehen werde. 

Indem wir die komplizierten Terme mit den trigonometrischen Funktionen durch die Teilfunktionen 
f 1 , f2 und f3 ersetzen, konnen wir die Schnittkurve kompakter darstellen 


One obtains exactly the same invtersection curve when solving the z-coordinate of the 
original parameterization Spi ralTorusX_u(b,m, t ,c|)) for ^ and then substitutes for cp. I don’t 
want to prove this here but I ask the interested reader to try this proof. 

The choice of this alternative representation has a second advantage: we can find appro- 
priate initial values for the Newton method which will be used later in my paper. 

Replacing the complicated expressions containing the trig functions by sub-functions f1, f2, 
and f3 it is possible to give the intersection curve a more compact form. 

#33: fl(b, t) :z 2 ■ SIN(2 -jt- t) ■ LN(b) ■ CLN(b) ■ SIN(2 -jt- t) + 2 -jt. C0S(2 -jt- t)) 

2 2 2 2 

#34: f2(b, m, t) := (m .COSCZ-T^.t) + D-LNCb) + 4.jr.m ■ SINC2 -jt- t) ■ (jt- SIN( 2 -fr. t) - LN(b) ■ C0S(2 -jt. t)) 

2 2 

#35: f3(b, m, t) :z (m + 1) ■ LN(b) ■ C0S(2 ■ jt- t) - Z-jr-m -SINCZ-JT-t) 
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#36: 


SCHNITTCb, m, t) = 


■ 2 t 

" flCb, 

t) 

m ■ b ■ 



f2(b, m, 

1 

. t) 2 


. 


m j 


t f3(b, m, t) 

2'LN(b)'m-b ■ 

f2(b, m, t) 


Schauen wir uns ein paar Schnittkurven fur Schneckenhauser an, die die von uns geforderte 
Bedingung erfiillen, namlich dass sie sich bei jeder Umdrehung selbst beriihren ohne sich dabei zu 
durchdringen. Die beiden Diagonalen stellen den Schnitt des Kegels dar, in dem sich die Geriistspirale 
windet. 

Let’s inspect some intersection curves which are satisfying the required condition by touching 
themselves at each rotation without interpenetration. The two lines represent the intersection 
of the cone which contains the winding base spiral. 


Example: 

SCHNITTC2 , 343263/1083221 , t) 

and 

-1/2 < t < 3/2 


Example: 

SCHNITTCe , 571877/1286241 , t) 

and 

-1/2 < t < 3/2 



Example: 

SCHNITTC3 , 7166353/14797132 , t) 

and 

-1/2 < t < 3/2 



Anhand der Grafiken erkennt man, dass beim Schneckenhaus vom Typ 2 die Bestimmung des passen- 
den m leider nicht so leicht ist wie beim Typ 1 . 

Beim Typ 1 liegen die Beriihrungspunkte zumindest auf der Y-Achse, so dass die Parameter (p der 
beiden Beriihrungspunkte analytisch ermittelt und in die Schnittkurven eingesetzt werden konnen. 
Dadurch reduziert sich dort das Problem auf eine einzige Gleichung. 

Hier sehen wir aber, dass beim Typ 2 die beiden Beriihrungspunkte weder auf der Y- bzw. X-Achse, 
noch auf den Kegelerzeugenden liegen. Die beiden Parameterwerte t der Beriihrungspunkte lassen sich 
hier leider nicht analytisch ermitteln und in die Schnittkurve einsetzen. Anstelle dessen miissen wir 
beim Typ 2 sowohl m als auch die beiden Werte fur t numerisch bestimmen. 
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Da nun die Werte von drei Variablen gesucht sind, miissen wir das Problem als ein System von drei 
Gleichungen formulieren deren Nullstellen wir bestimmen miissen. 


Wie sehen diese Gleichungen aus? Zum einen miissen natiirlich die beiden Beriihrungspunkte die- 
selben Koordinaten haben. Die Variablen r und s bezeichnen bier die beiden verschiedenen Parameter 
t der Beriihrungspunkte. Daraus ergeben sich die ersten beiden Gleichungen. 

Fiir die X- und Y-Koordinate der Beriihrungspunkte erhalten wir - mit r> s\ 

Inspecting the graphs we can see that finding the appropriate m for type 2 shells is not as 
easy as it was for the type 1 shells. 

The tangency points for type 1 are at least lying on the y-axis which enables finding the 
parameter q) of both points analytically and then substituting in the intersection curves. The 
problem is reduced to one single equation. 

But here for type 2 both tangency points are neither lying on one of the two axes nor on the 
surface line of the cone. Both values for t of the tangency points cannot be determined 
analytically for substitution into the intersection curve. The values for m and for both t can 
only be found numerically. 

As values for three variables are requested, the problem must be defined as a system of 
three equations in order to determine three zeros. 

How do these equations look like? First of all both tangency points must have the same 
coordinates. Variables rand s indicate the two different parameters t. This leads to the first 
two equations. 

For the x- and y-coordinate of the tangency points with r> s: 


#37; 


#38; 


0 = (SCHNITTCb, m, r)) - (SCHNITT(b, m, s)) 

2 2 

0 = (SCHNITTCb, m, r)) - (SCHNITT(b, m, s)) 

1 1 


Gleichzeitig soil sich die Schnittkurve an den beiden Benihrungspunkten nicht selbst schneiden. Als 
Konsequenz daraus stehen die Tangenten der Schnittkurve an den beiden Beriihrungspunkten parallel 
zu einander. Auch das lasst sich als eine Gleichung formulieren: 

Simultaneously the intersection curve shall not intersect itself in both tangency points. We 
will denote this in an equation expressing that the tangents in these points must be parallel: 



fd ^ 

d 

fd 'l 

#40: 0 = 

— (SCHNITTCb, in, r)) 

(SCHNITTCb, m, s)) - 

— (SCHNITTCb, m, r)) 


^dr 1, 

ds 2 

^dr 2, 


(SCHNITTCb, tn, s)) 

1 

Damit haben wir ein Gleichungssystem formuliert. Gesucht ist ein Tripel m, r und 5 fiir das gilt: 
So we have set up a system of equations for unknown m, r, and s as follows: 


a 


2 r 
m ■ b 


fl(b, r) 1 ^ 

+ 1 

f2(b, m, r) 2 

.. m j 


2 s 
m ■ b 


f fl(b, s) 1 

+ 1 

fZ(b, m, s) 2 

^ m j 


#41: 


p 40 


Piotr Trebisz: A Mathematical Model for Snail Shells 


D-N-L#82 


r f3(b, m, r) s f3Cb, m, s) 

#42: 0 = 2.LNCb).m-b 2.LN(b).m.b 


#43: 


1 


(dr 

[■ ‘ 1 


2.LN(b).m.b 


f2(b, m, r) 

fl(b. r) 1 

+ — 1 

f2(b, m, r) 2 

m 


r f3(b, m, r) 



' d " 


ds ^ 

JJ. 

* 


f2(b, m, s) 

s f3(b, m, s) 

2' LN(b) ‘m* b 


f2(b, m, s) 


]-& 


f2(b 


, m, r) y 

. m. r) 


11 

d 

"2s" 


. — 

m «b • 

JJ 

ds 






fl(b, s) 1 

+ — 1 

f2(b, m, s) 2 

m 


Die drei Gleichungen konnen waiter vereinfacht warden zu: 

These equations can be “simplified” further to: 


r 

0 = b - 

#47: 

^ f iCb , 

r) 

1 ^ 
+ - 1 

s 

- b ■ 

H 

M- 



S> 

1 ^ 
_ 1 

fZ(b, m 

, r) 

Z 

m j 

fZ(b, m 

, s) 

z 

m j 

r 

#48: 0 = b ■ 

f3(b, m, 

r) 

s f3Cb, 
■ 

m, s) 




fZ(b, m, 

r) 

f2(.b, 

m, s) 





d f r f fl(b , r) 


dr 


#49: 


b - 




+ 1 

fZ(b, m, r) 2 

K m JJ 


d f s f3(b, m, s) 
b - 
ds 


fZCb, m, s 


7] 


0 = 


d 

dr 



d 

" s 

" flCb, s) 

1 

f r f3(b, m, r) ^ 

— 

b - 


+ 1 

b 

ds 


fZCb, m, s) 

Z 

L fZ(b, m, r) J 




m jj 


r 

b 


s 

b 


Dieses Gleichungssystem ist hochst nichtlinear. Deshalb rniissen wir uns, genauso wie beim 
Schneckenhaus vom Typ 1, des Newton-Verfahrens bedienen, um dieses Gleichungssystem numerisch 
zu losen. Gesucht sind die Nullstellen der folgenden Funktionen: 

This system is not linear at all. So we have to apply - like we did for the type 1 shells - the 
Newton-Raphson method in order to obtain the zeros of the following functions numerically: 


]l(b, m, r, s) 

#50: 

r 

:~ b - 

^ flCb , 

r) 

1 "1 

+ - 1 

s 

_ . 

" f iCb , 

s) 11 

I - 1 


fZ(b, m, 

, r> 

z 

m j 


fZCb, ni: 

, s) Z 

rti j 

#51: ]Z(b, m, r, s) 

r 

— 

f3(b, m, 

r> 

s f3(b, 
— b ■ 

m, s) 




f2Cb, m, 

r> 

f7Cb, 

m, s> 




# 52 : 


]3(b, m, 


s) := 


d 

dr 


f r ( flCb, r) 
b ■ 


fZ(b, m, r> 




r 

b 


1 

Z 

m 




1 


JJ 


d 

ds 


s f3(b, m, s) 

b 

fZ(b, m, s} J 


b 




d 

r s 

^ flCb, s) 

1 

d 

f r f3(b, m, r) 

— 

b - 


^ 1 

— 

b . 

ds 


fZ(b, HI, s) 

Z 

dr 

L fZ(b, m, r) J 




m Jj 


r 


s 


b 


b 
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Dies ist eine Iteration des mehrdimensionalen Newton-Verfahrens. 

This is an iteration of the multi-dimensional Newton method. 

d d 


#53: [ni2, T2, S 2 ] - [mi, Pi, Si] 


— 

31(b, 

mi. 

Tl, Si) 

— 31(b, 

mi. 

Tl, 

Sl) 

— 31(b, 

mi. 

Tl, Sl) 

dm 




dp 




ds 



d 




d 




d 



— 

J2(b, 

mi. 

Ti, Si) 

— 32(b, 

mi. 

r 1 , 

Si) 

— 32(b, 

mi. 

Tl, Sl) 

dm 




dp 




ds 



d 




d 




d 



— 

J3(b, 

nil. 

Tl, Si) 

— 33(b, 

mi. 

ri. 

Sl) 

— 33(b, 

mi. 

Tl, Sl) 


dm 


dr 


ds 


[Jl(b, mi, Pi, Si), J2(b, mi, Pi, Si), J3(b, mi, Pi, Si)] 


Eine aquivalente, aber effizientere und kompaktere Darstellung fiir die Iteration lasst sich mit einer 
Jacobi-Matrix erzielen. 

An equivalent but more efficient and more compact form of this method can be achieved 
using a Jacobi matrix. 


#54: 3AC0BI_MATRIX(b, mi, Pi, Si) := 

nil, Ti, Si) 

nil, Ti, Si) 
mi. Pi, Si) 

#55: [ni2, P 2 , S 2 ] = [mi, Pi, Si] - f 


d d 

— 31(b, nil, Ti, Si) — 31(b, 

dm dp 

d d 

— J2(b, mi, Pi, Si) — 32(b, 

dm dp 

d d 

— J3(b, nil, Ti, Si) — 33(b, 

dm dp 


LREDUCE(3AC0BI_MATRIX(b, mi. 


d 


nil. 

Tl, 

Sl) 

ds 

31(b, 

nil. 

Tl, Sl) 

31(b, 

nil. 

Tl, 

Sl) 

d 

ds 

32(b, 

nil. 

Tl, Sl) 

32(b, 

nil. 

Tl, 

Sl) 

d 

ds 

33(b, 

nil. 

Tl, Sl) 

33(b, 


, Si)) 


Die drei Funktionen Jl, J2 und J3 sowie ihre Ableitungen nach m, r und a sind sehr komplex, deshalb 
rniissen wir sie in handliche Teilfunktionen zerlegen um ihre Komplexitat zu bandigen, ganz nach dem 
Prinzip "Teile und Herrsche". 

The three functions J1 , J2, and J3 and their derivatives wrt m, r, and s are very complex. 
Thus we have to decompose them into comfortable sub-functions following the principle 
“Divide et impera” = “Divide and Conquer”. 

#56: [vl :€ Vector, v2 :€ Vector, v3 :c Vector] 

#57: [Fr :g Vector, Fs :€ Vector, Gr :€ Vector, Gs :€ Vector] 

Als erstes definieren wir uns eine Funktion F, die aus den Teilfunktionen fl, f2 und f3 und ihren 
Ableitungen nach m und t besteht. 

Firstly we define a function F consisting of the sub-functions f1 , f2, and f3 and their 
derivatives wrt m and t. 
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#58: FCb, m, t) := 


flCb, t) 


d 

— flCb, t) 
dm 


d 

— flCb, t) 
dt 


d d 

fl(b, t) 

dm dt 


d 

fZ(b, m, t) — f2(b, m, t) 
dm 


d d d 

— f2(b, m, t) f2(b, m, t) 

dt dm dt 


f3(b, m, t) 


d d 

— f3(b, m, t) — f3(b, m, t) 
dm dt 


d d 

f3(b , m, t) 

dm dt 



flCb, t) 
f2(b, m, t) 
f3(b, m, t) 


Nach dem Baukastenprinzip konnen wir daraus die Teilkomponenten der Jacobi-Matrix zusammen 
setzen. Die Vektoren vl, v2 und v3 reprasentieren bier die Teilfu nk tionen fl, f2 und f3 und ihre 
Ableitungen nach m und t. 

Ich babe die Terme nach dem Homer-Schema definiert. Dadurch ist die Darstellung kompakter und 
die Auswertung effizienter. Ausserdem wirkt sich das Homer-Schema positiv auf das Verhalten von 
Rundungsfehlem aus. 

According to the modular system we can assemble the subcomponents of the Jacobi matrix. 
Vectors v1, v2 and v3 are representing the sub-functions f1 , f2 and f3 and their derivatives 
wrt m and t. 

I defined the expressions following the Horner procedure. This makes presentation more 
compact and evaluation more efficient. Additionally this treatment has a positive effect on the 
behaviour of rounding errors. 

Expressions #60 and #61 are so bulky that I don’t print them here. I’d like to invite you 
inspecting the original DERIVE file. 

This are G(b,m,vl,v2,v3) and G(b,m, (F(b,m,t))i,(F(b,m,t)) 2 ,(F(b,m,t) 3 )). 

You can see only the first part of the first module. Josef 


#60: G(b, m, vl, v2, v3) 


vl vl ■ v2 

11 12 2 

v2 2 2 3 

1 m v2 m 


v3 


v3 


v2 


v2 -v3 
2 1 


v2 


v2 


vl + vl ■ 
3 1 


v2 

3 


LN(b) 


v2 


1 


f 1 


■LNCb) 


v2 


v3 + v3 ■ 
3 1 


v2 


LN(b) - 


v2 


v2 


1 J 


Leider ist DERIVE damit iiberfordert diesen Ausdmck auszuwerten weil der Speicher nicht ausreicht. 
Doch mit etwas Ubung und durch Anwendung gmndlegender Ableitungsregeln konnen die interes- 
sierten Leser die Richtigkeit meiner Zerlegung auch per Hand verifizieren. 

Damit sind wir nun endlich in der Lage aus den einzelnen Komponenten die Jacobi-Matrix 
zusammenzusetzen. 

Unfortunately is DERIVE unable to evaluate this expression - lack on memory. But with 
some routine interested readers can verify my decomposition manually applying basic 
differentiation rules. 

Finally we are able to compose the Jacobi-matrix of its components. 
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#62; 

[Fr := F(b, m, r), Fs : 

F(b, m, s)] 








f Gr := G(b , m, Fr , Fr , 

Fr ), Gs := 

G(b, m, F 

■s , Fs , Fs 

)l 




#63 : 

L 1 

1 

3 


12 3 

J 







r 

s 




r 




b -Gr 

- b ■Gs 




b ■Gr 




1.2 

1,2 




1,3 




r 

s 




r 

#64: 

DACOBI_MATRIX(b , m, r, s) :r 


b -Gr 

- b ■Gs 




b ■Gr 




2,2 

2,2 




2,3 



Gr -Gs + Gr -Gs 

- Gr ■Gs 

Gr 

■Gs 

Gr ■Gs 

; - Gr ■Gs 



L 1,4 2,3 

1,3 2,4 

2,4 1,3 

2, 

3 1, 

4 1,5 

2,3 2,5 1,3 


s 

r s 








- b -Gs 

b -Gr - b 

■Gs 







1.3 

1,1 

1,1 







s 

r s 








- b -Gs 

b -Gr - b 

■Gs 







2,3 

2,1 

2.1 







Gr -Gs - Gr -Gs 

Gr -Gs - Gr 

■Gs 







1,3 2,5 2,3 1,5 

1,3 2,3 

2,3 1,3 . 








■ d 

d 


d 






— 31(b, m, 

r. s) — 

31(b, m, r, s) 

— 

31(b, 

m, r, s) 

31(b, m, r, s) 



dm 

dr 


ds 






d 

d 


d 




#65: 

3AC0BI_MATRIX(b , m, r, s) 

= 

— 32(b, m, 

r, s) — 

32(b, m, r, s) 

— 

32(b, 

m, r, s) 

32(b, m, r, s) 



dm 

dr 


ds 






d 

d 


d 






— 33(b, m. 

r, s) — 

33(b, m, r, s) 

— 

33(b, 

m, r, s) 

33(b, m, r, s) 



. dm 

dr 


ds 





Unser Newton-Verfahren ist fast fertig. Was uns jetzt noch fehlt sind die Startwerte fur m, r und 5 und 
die Einschrankung von r und 5 auf ein Losungs-Interval. 

Damit das Newton-Verfahren zu einem sinnvollen Ergebnis konvergiert sind verniinfitige Startwerte 
fur die Variablen unerlasslich, die nicht allzuweit entfemt sind von den gesuchten Werten. Doch 
welche Startwerte fur m, r und 5 sollen wir wahlen? Diesbeziiglich erweist sich die alternative 
Parametrisier-ung Schnecke_U (b , m , t , c|)) einmal mehr als sehr hilfreich. Wenn wir uns den Aufbau 
dieser Funktion genauer anschauen, dann stellen wir fest, dass sie viele Gemeinsamkeiten mit dem 
Schneckenhaus vom Typ 1 besitzt. Unser raumliches Vorstellungsvermogen ist jetzt gefragt! 

Typ 1 basiert auf einer Spirale in der YZ-Ebene. Die lokale mitdrehende X-Achse ist dort ein Vektor 
der senkrecht zur YZ-Ebene steht und der in Richtung der positiven globalen X-Achse zeigt. Die 
lokale mitdrehende Y-Achse steht senkrecht zur Normalen der YZ-Ebene und zeigt dabei weg vom 
Zentrum der Spirale. 

Ganz ahnlich verhalt es sich bei der Funktion Schnecke_L) (b , m , t , (()) . Sie basiert auf einer Spirale die 
sich auf der Oberflache eines Kegels windet. Ihre lokale mitdrehende X-Achse steht auch senkrecht 
auf der Oberflache des Kegels. Das lasst sich leicht zeigen. 

The Newton-procedure is almost ready. What we are missing are the initial values for m, r, 
and s and the restriction on a solution interval for rand s. 

In order to obtain a reasonable result from the Newton method it is a must to have appro- 
priate initial values for the unknown which are not too far from the requested values. But 
which initial values should we choose? Once more the alternative parameterization 
Schnecke_U(b,m, t ,(()) proves to be very helpful. Inspecting the structure of this function we 
can recognize many similarities with the type 1 snail shells. This is a challenge for our spatial 
sense! 
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Type 1 is based on a spiral in the yz-plane. The local co-rotating x-axis is a vector perpend- 
icular to the yz-plane directed to the positive global x-axis. The local co-rotating y-axis is 
perpendicular to the normal of the yz-plane directed out of the spiral centre. 

Very similar is the situation with the function Schnecke_U(b,m,t,(t)). It is based on a spiral 
winding on the surface of a cone. Its local co-rotating x-axis is perpendicular to the surface of 
the cone. This can easily be shown. 

z 

#66: LINIECb, m, t, ra -t + t> • SNAILfb, m, t) = 0 


2 d 

#67: LINIECb, m, t, m -t + t) • — SNAILCb, tti, t) = 0 

dt 


Die lokale mitdrehende X-Achse zeigt in Richtung der positiven globalen X-Achse 


The local co-rotating x-axis shows in the direction of the positive global x-axis. 

2 

#68: (LINIECb, m, t, m -t t)) > 0 

1 

Die lokale mitdrehende Y-Achse steht senkrecht zur Normalen der Kegeloberflache und zeigt weg 
vom Zentrum der Spirale 


The local co-rotating y-axis is perpendicular to the normal of the surface of the cone and is 
directed out of the centre of the spiral. 


#69: 


Scbnecke_U 

b, m, t, 

> 

Schnecke_U 

fb, t, - -11 


L 2 J 



L 2 )\ 


Im Grunde ist Schnecke_U(b,m, t ,c|)) also nichts anderes als ein Schneckenhaus vom Typ 1, dessen 
YZ-Ebene zu einem Kegel zusammengefaltet wird. Es ist daher nahe liegend fiir m, r und a dieselben 
Startwerte zu verwenden wie bei der Berechnung von m fiir das Schneckenhaus vom Typ 1 . 

Basically Schnecke_U(b,m, t ,((>) is nothing else than a type 1 snail house, with its yz-plane 
folded to a cone. Hence it seems to make sense using the same initial values for m, rand s 
which we had calculating the m for the type 1 snail house. 

M_YZ(b,n) is defined like for Trebisz Spiral Type 1! (You can find it in the file!) 


#71: mt> := M_YZ(b , n) 


#72: 


= 2-ATAN 


r 2 2 

(fiio -F l)-LNCb)- C0TC2-TT- To) - 2-TT-mQ 

2 2 2 2 
V ^CCtiio + l)~LN(b) 4'7 t mo ) 


#73 : -r = 2 -ATAN 

2 


r 2 2 

(mo + l)-LNCb)"COT(2-7T- So) - 2~n mo 

2 2 2 2 
/((mo + l) LN(b) + 4'7t nio ) 


Fiir "ro" und "sq" bedeutet das: / This is for "ro" ^nd "Sq": 


# 74 : 




ATAN 


2 

^ VCCmo 


2 ^ 
Cttio + l)-LN(b> 

2 2 2 2 
+ l)-LN(b) + 4-tt "Hio ) - 2-jr-riio J 


2-Ji 


Tq - 1 - 
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ATAN 


#75: 


(ttio + l}-LNCb) 


5f> - 


Z Z Z 2 Z 

yCCnio + l)"LNCb) + 4-tt -mo ) + Z-rr-mo 


Z-7T 


Mit der 1 vor dem negativen Arcus-Tangens-Term in der Funktion fiir ro wird dem Umstand 
Rechnung getragen, dass ro ungefahr eine Umdrehung weiter liegt als sq . Die Funktion fur ro ist 
allerdings noch nicht ganz korrekt. Der Arcus-Tangens-Term von ro ist namlich nicht stetig. 

Wenn der Nenner in dem Arcus-Tangens-Term . . . 

The 1 in front of the negative arctan-expression is due to the fact that ro is approximatively 
one rotation in front of Sq. The function for "ro" is not yet completely correct. The arctan- 
expression is not continuous. 

When the denominator in this expression ... 

z z z z z 

#76: JCCmo l)-LNCb) + 4'TT -mo ) - Z-yr-mo 

. . . negativ wird, dann wird ro sprunghaft um 1/2 groBer. Um das zu kompensieren fiihren wir einen 
Korrekturterm ein. 

... becomes negative, then ro increases with a jump by 1/2. In order to compensate this we 
introduce a correcting expression. 

(2 ^ 
(mo + l)-LN(b) 

ATAN 

#77: 2 222 2 2 222 2 

3 + SIGNCACnio + l)-LN(b) + 4 -tt -mo ) - 2-7r-mo ) L /((mo + l)-LN(b) + 4 -tt -mo ) - 2-7r-mo > 

To = 

4 2-tt 


Es empfiehlt sich unsere theoretischen Annahmen iiber die Startwerte grafisch zu iiberprufen 


It is recommended to double check our theoretical assumptions about the initial values 
graphically. 


Example: 

SCHNITT(2,mo,t) a 
SCHNITT(2,mo,So) a 
SCHNITT( 2,mo, To) 

with -1/2 <t< 3/2 


Example: 

SCHNITT(e,mo,t) a 
SCHNITT(e,mo,So) a 
SCHNITT( e,ttio, To) 




with -1/2 <t< 3/2 
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Example: 

SCHNITT(3,mo,t) a 
SCHNITT(3,ttio,So) a 
SCHNITT( 3,ttio, To) 

with -1/2 <f<3/2 



Wie man sieht, liegen die Startwerte in der Tat nah bei den tatsachlichen Werten. Was jetzt noch fehlt 
ist die Einschrankung von r und s auf ein Losungsinterval. 

Dies ist notwendig da die Schnittkurve eine durchgehende Linie ist. Es besteht die Gefahr dass r und s 
entlang der Schnittkurve wandem bis sie sich irgendwo in der Mitte treffen. Wir miissen deshalb 
sicherstellen, dass r und s in den beiden Schlaufen verbleiben wo sich die beiden Beriihrungspunkte 
befmden. Die beiden Schlaufen beginnen und enden an den Nulldurchgangen der Y-Koordinate der 
Schnittkurve. 

We can see that the initial values really are very close to the actual values. What is missing 
is the restriction of rand s on a solution interval. 

This necessary because the intersection curve is one continuous line. There is the danger 
that rand s are moving along the intersection curve until they are meeting anywhere in the 
middle. So we have to make sure that rand s will remain in those both loops where the 
tangency points are lying. Both loops start and end at the zeros of the y-coordinate of the 
intersection curve. 


0 = (SCHNITTCb, m, t)) 
#78; 2 


#79; 0 = 2-LN(b)-m-b ■- 


2 2 
(m -I- l)-LN(b)-COS(2-7T- 1) - 2-7r-m •SIN(2-7T-t) 


2 2 2 2 

(m •COS(2-7T*t) + l).LN(b) + 4-7T-m -SIN(2*7r- 1). (7r-SIN(2-7T-t) - LN(b).COS(2-7T- 1)) 


Eine Reihe moglicher Losungen sind: 

Some possible solutions are: 


ATAN 


#80: 


(m -b l).LN(b) 


2 -TT-ni 


ATAN 


(m -b l).LN(b) 


2 -TT-m 


ATAN 


(ni -b l).LN(b) 


2 -n-m 


2 -TT 


ATAN 


2 -TT 


2 -TT 


(m -b l).LN(b) 


2 -TT -m 


-b V t4 

2 


-b 1 


2 -TT 


Durch Ausprobieren konnen wir herausfinden in welchen Grenzen die beiden Schlaufen liegen. 

By experimenting we can find out the boundaries for the two loops. 


D 
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Example: 


b = e A m = 571877/1286241 a -1/2 < t < 3/2 and 
SCHNITTCb, m, ti) -> Rot / Red 

SCHNITTCb, m, tz) -> Gelb / Yellow 

SCHNITTCb, m, ts) -> Griin / Green 

SCHNITTCb, m, 14 ) -> Blau / Blue 



Now we know the restrictions for the solution intervalls for r and s. 


A 2 

Cm + l)-LNCb) 


ATAN 


#81: 


Z-TT-m 


r 2 A 

(m -If l)-LNCb) 


ATAN 


Z-TT-m 


#8Z: 


Z-TT 

r 2 ^ 

(m -I- l)-LNCb) 

Z 

Z-TT-m 


+ < r < 

Z 


ATAN 


Z-TT 

r 2 ^ 

Cm -If l)-LNCb) 

Z 

Z-TT-m 


1 


ATAN 


Z-TT 


Z-TT 


Und was passiert wenn im Laufe des Newton- Algorithmus r oder a das Intervall verlassen? Ganz 
einfach, wir setzen es auf einen neuen Startwert mit dem aktuellen m zuriick. 

Mit all den Erkenntnissen, die wir gewonnen haben sind wir nun endlich in der Lage den Newton- 
Algorithmus, der uns das passende m und die beiden Umdrehungszahlen r und a der beiden Beriih- 
rungspunkte berechnet, in ein DERIVE-Programm zu gieBen. Das Programm berechnet m, r und a auf 
n Stellen genau 

And what will happen when r or s are leaving the solution interval during the Newton- 
algorithm? Easy answer: we start again with the actual m as new inital value. 

Applying all our freshly acqired insights we are now able to bring the Newton algorithm which 
shall calculate the appropriate m and the two parameters rand s of both tangency points in 
the form of a DERIVE program. The following program returns the values of m, r, and s with 
an accuracy of n decimal places. 

In his last mail Piotr gave a smaller interval for rand s together with the respective explication. 

You can find this improvement in the dfw-file. 
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#8J: I(b, m) := 


ATAN 


(m + l)'LN(b) 


2 'JT' m 


ATAN 


ATAN 


2 ■ JT 


Cm + l)'LN(b) 


2 'IT' m 


2 -IT 


(m + l)'LNCb) 


2 -JT' m 


2 -jr 

f 2 

(m + l)'LN(b) 


+ 1 


ATAN 


2 ■JT' m 


2 -JT 


ATAN 


2 2 2 2 2 
3 + SIGNC-/C(m + l)'LN(b) + 4 -jt -m ) - 2-iT'm ) 


(m + l)'LNCb) 


2 2 


^ ^i^CCrn + l)'LN(b) + 4 -?t -m ) - 2'7T'm j 


2 'JT 


ATAN 


(m + l)'LN(b) 


2 2 2 2 2 
7C(m + l)'LN(b) + 4-]nr -m ) + 2'JT-m j 


2 ■?? 


H(b, m, r, s, n) := 

Prog 

[Br :g Real (1, b), Bs :e Real (1, b), RS :e Vector] 

[Fr :e Vector, Fs :e Vector, Gr ;e Vector] 

[Gs ;e Vector, 3M ;e Vector, AM ;e Vector] 

WRITE CAPPENDC"m APPROX(STRINGCm) , n))) 

RS := APPROXCiCb, m), 2-n) 

If RSil^l < r < RS4,1|2 
r := RSil|3 

If RSx2il < s < RS4.2 x 2 
s := RS| 24,3 

[Br APPROXCb^r, 2-n), Bs APPROX(b^s, 2-n)] 

[Fr := APPROXCFCb, m, r), 2-n), Fs := APPROX(FCb, m, s), 2-n)] 

#S4; Gr := APPROX(GCb, m, Fr4,l, Fr|2, Fr|3), 2-n) 

Gs APPROXCGCb, m, Fs|l, Fs|2, Fs|3), 2-n) 

3M:=[?, ?, ?, ?; ?, ?, ?, ?; ?, ?, ?, ?] 

[IMil^l := Br-Gr4,l|2 - Bs'Gs|l4,2, 3M|1|2 := Br-Gr4,l4,3] 

[3Mj.li3 := - Bs'Gsxlx3, 3Milx4 := Br-Grilil - Bs-Gsxlil] 

[3Mx2il := Br-Gri2x2 - Bs-Gsi2x2, 3Mx2x2 := Br-Grx2x3] 

[3Mx2x3 := - Bs'Gsx.2x3, 3Mx2|4 := Br'Grx.2xl - Bs-Gsx2xl] 

3Mx3il := Grxlx4- Gsx2x3 + Grxlx3 ■ Gsi2x4 - Grx2x4- Gsxlx^ - Gr x2x3 ■ Gsxlx4 
3Mx3x2 := Grxlx^ ■ Gsx2x3 - Grx2x5 ■ Gsxlx^ 

3Mx3x3 := Grxlx^ ■ Gsx2x3 - Grx2x3 ■ Gsxlx^ 

3Mx3x4 := Grxlx^ ■ Gsx2x3 - Grx2x3 ■ Gsxlx^ 

AM := [m, r, s] - (R0W_REDUCEOM))xx4 
RETURN APPROXCAM, n) 

^LXCb, n) := 

Prog 

[mo :e Rea] (0, oa), rg :e Real (0, 1), sg :e Rea] (0, 1)] 

V :e Vector 

#S5: mg - M_¥ZCb, n) 

rg := APPROXCaCb, mg))xlx^, n) 
sg := APPROXCCKb, mg))x2x3, n) 

V := ITERATES (H (b , m, r, s, n), [m, r, s], [mg, rg, sg]) 

RETURN VxDIM(V) 


Damit sind unsere konischen Schneckenhauser vom Typ 2 endlich fertig! Man muss nur das passende 
m in die entsprechenden Formeln einsetzen; 

Now our conic snail shells are finally ready. We only have to substitute the appropriate value 
for m into the respective formulae: 


Spi ralFlacbeX_L(b , n, t, cjs) := Spi ralToriJsX_LCb , (M_XCb, n)) , t, tp) 
#S6; 1 

Spi ralFlacbeX_UCb , n, t, 0) := Spi ra]Tor JsX_U(b , CM_X(b, n)) , t, tp) 

1 
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Die Differenz zwischen r und s entspricht genau einer Windung, also der Umdrehungszahl nach der 
sich der Mantelkreis des Schneckenhauses selbst beriihrt. Damit konnen wir eine dritte Variante des 
Schneckenhauses definieren, das nun nach der Windungszahl parametrisiert ist. 

The difference between rand s corresponds to exact one twist, i.e. the winding number after 
which the surface of the snail house touches itself. So we are able to define a third form of 
the snail house which is parameterized according to ist twist number. 


#SS; LI i£ Vector 

Spi ralTorusX_W(b , u, t, up) := Spi ralTorusX_U(b , - u ), up) 

#S9: 12 3 

#90; Spi ralFlacheX_W(b , n, t, up) := Spi ralTorusX_W(b , M_X(b, n), t, tp) 

Unsere Anstrengungen werden mit asthetischen 3D-Grafiken der Schneckenhauser belohnt. 

Our efforts are rewarded with esthetic 3D plots of the snail houses. 


Base 2: Spiral FI acheX_w(2, 3) (-3 s s ^ 3, 100 panels; 0 ^ f ^ 27i, 50 panels) 



Base 6: Spi ral FlacheX_W(6, 3) (-2 < s < 0, 100 panels; 0 < f < 27i, 50 panels) 


y 
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Basis e mit genau einer Windung: 
Base e with exact one twist: 

Spi ra1 F1acheX_W(e, 3) 

(-1 < s < 30, 50 panels; 

0 < t < 271, 50 panels) 



Important note: As mentioned earlier in DNL#81 only the mth-file can be used to plot the 
surfaces. \Ne - Piotr and I - don’t have any idea why the functions don’t work when stored in 
dfw-format. I wrote to Albert Rich and hopefully I can deliver an explanation next time. Josef 


I couldn’t resist creating my own snail house. I am very curious, so let’s have a look into the 
interior of the home of a snail! 

Spiral FI acheX_W(l. 25, 3) (-3 S s ^ 3, 120 panels; 0 ^ t ^ 7i, 30 panels) 

(Box dimensions: 0,3; -1 .5,1.5; -1 .5,1.5) 





